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Abstract. From the correlation structure of transverse momentum p t in relativistic nuclear 
collisions we observe for the first time temperature/velocity structure resulting from low- 
Q 2 partons. Our novel analysis technique does not invoke an a priori jet hypothesis. p t 
autocorrelations derived from the scale dependence of (pt) fluctuations reveal a complex parton 
dissipation process in RH1C heavy ion collisions. We also observe structure which may result 
from collective bulk-medium recoil in response to parton stopping. 
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Q^' 1. Introduction 

Central Au-Au collisions at RHIC may generate a color-deconfined medium (quark-gluon plasma 
or QGP) pQ. Some theoretical descriptions predict abundant low-Q 2 gluon production in the 
early stages of high-energy nuclear collisions, with rapid parton thermalization as the source 
of the colored medium [2 El 13- Nonstatistical fluctuations of event- wise mean pt {pt) El 
may isolate fragments from low-Q 2 partons and determine the properties of the corresponding 
medium. A recent measurement of excess {pt) fluctuations in Au-Au collisions at 130 GeV 
revealed a large excess of fluctuations compared to independent-particle pt production . 

In this paper we describe the event-wise structure of transverse momentum pt produced in 
relativistic nuclear collisions at RHIC. We discuss the role of low-Q 2 partons as Brownian probe 
particles in heavy ion collisions. We compare joint autocorrelations on (77, 0) to conventional 
leading-particle techniques for parton fragment analysis. We present experimental evidence 
from mean-pt fluctuations and corresponding pt autocorrelations for local temperature/velocity 
structure in A-A collisions which can be interpreted in terms of parton dissipation in the A- 
A medium and same-side recoil response of the bulk medium to parton stopping. Finally, we 
review the energy dependence of mean-p^ fluctuations from SPS to RHIC and its implications. 



2. Low-Q 2 partons as Brownian probes 

In 1905 the microscopic structure of ordinary matter was addressed theoretically by Einstein, 
who introduced the concept of a (Brownian) probe particle large enough to be observed 
visually, yet small enough that its motion in response to the molecular dynamics of a fluid 
might also be observed [Zj. Those two constraints specified the one- micron probe particles 
used by Jean Perrin to confirm molecular motion in fluids [HI Ej- The Langevin equation 
a stoch(t) + Omcs(^) models the motion of a Brownian probe in a thermalized fluid 



medium of point masses qualitatively smaller than the probe particle jlUUllj . The accelerations 
are gaussian-random with zero mean; a stoc h(t) is isotropic and a mcs (t) _L v(t) (and oc v). The 
first term models collective dissipation of probe motion (viscosity), the second models individual 
probe collisions with medium particles and the third simulates multiple Coulomb scattering of a 
fast probe particle. A solution of that equation for unit initial speed in the x direction starting 
at the (x, y) origin is shown in the first two panels of Fig. ^ Speed is dissipated with time, 
leading to equilibration with the medium: fluctuations of velocity about zero and random walk 
of the probe. An example of such motion is shown in the third panel: an electron track in a 
time projection chamber exhibits multiple Coulomb scattering along its trajectory, terminating 
in random walk represented by the circled ball of charge at the trajectory endpoint |12j . 
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Figure 1. Left panels: Solution to the Langevin equation for a point mass with unit initial 
velocity in the positive x direction illustrating thermalization. Right panel: two views of an 
electron track in a TPC, with charge ball from random walk at track endpoint (circled) |12j . 

In 2005 we seek the microscopic structure and local dynamics of the QCD medium formed 
in RHIC heavy ion collisions. The point-mass concept of Einstein's Brownian probe must be 
extended to partonic probes, possibly with internal degrees of freedom and experiencing complex 
non-point interactions with medium degrees of freedom. This problem requires novel analysis 
techniques closely coupled to the Langevin equation and its associated numerical methods. The 
analog in heavy ion collisions to Einstein's Brownian probe is the low-Q 2 parton, visualized for 
the first time by methods presented in this paper. In contrast to Einstein's notion of a particle of 
exceptional size observed indefinitely in equilibrium with microscopic motions of a thermalized 
particulate medium, the QCD Brownian probe is identical to medium particles but possesses 
an exceptional initial velocity relative to the medium with which it interacts for a brief interval: 
do probe manifestations in the hadronic system reveal 'microscopic' degrees of freedom of the 
medium, is the medium locally or globally equilibrated, what are its fluid properties? 



3. Joint autocorrelations vs conditional distributions 

Conventional study of QCD jets in elementary collisions is inherently model-dependent. 
Scattered partons with large transverse momentum are associated individually with 
concentrations of transverse momentum or energy localized on angle variables (??,</>). In 
heavy ion collisions, where such identification is impractical, jet studies are based on a high- 
Pt 'leading particle' which may estimate a parton momentum direction and some fraction of its 
magnitude. The leading-particle momentum is the basis for two-particle conditional distributions 
on transverse momentum and angles. Those distributions reveal medium modifications to parton 
production and fragmentation as changes in the single-particle pt spectrum (Raa) and in the 
fragment-pair relative azimuth distribution (away-side jet disappearance), referred to collectively 
as jet quenching [T3"j . The leading-particle approach is based on perturbative concepts of parton 
hard scattering as a point-like binary interaction and parton energy loss as gluon bremsstrahlung. 
We can then ask how the medium is modified by parton energy loss and what happens to low-Q 2 



partons, in a Q 2 regime where the pQCD assumption of point-like interactions breaks down, 
where the parton may have an effective internal structure. In other words, how can we describe 
parton dissipation as a transport process, including bulk-medium degrees of freedom? 

To access low-Q 2 partons we have developed an alternative analysis method for jet correlations 
employing autocorrelation distributions which do not require a leading- or trigger-particle 
concept. The autocorrelation principle is illustrated in Fig. |2J Projections of the two- 
particle momentum space of 130 GeV Au-Au collisions onto subspaces (771,772) and (0i,02) (left 
panels) indicate that correlations on those spaces are approximately invariant on sum variables 
7?s = 7?i+7?2 and <£>■£ = 0i + 02 , in which case autocorrelations on difference variables 77 a = 771—772 
and a = 01 — 02 retain nearly all the information in the unprojected distribution |14j . 
The autocorrelation concept was first introduced to solve the Langevin equation, to extract 
deterministic information from stochastic trajectories. In time-series analysis the autocorrelation 
of time series f(t) is A(t) = ^ (*) + r ) ^> wnere difference variable r = t\ — ti is the 

lag. For a stationary distribution ( f(t) correlations statistically independent of absolute time) 
the autocorrelation represents a projection by averaging of all the information in f{t). 




Figure 2. Left panels: Two-particle distributions on 77 and from central Au-Au collisions at 
130 GeV illustrating stationarity along the sum axis. Right panels: Joint autocorrelations on 
(t?a,<M) from p-p collisions (left) and Au-Au collisions (right). 

The same principle can be applied to ensemble-averaged two-particle momentum distributions 
which are approximately invariant on their sum variables [El- Distributions on angle space 
(7/1, 772, 01, ^2) can be reduced to joint autocorrelations on difference variables (t7a,0a)- For 
example, joint autocorrelations in the right-most two panels of Fig. |2] correspond to (771,772) and 
(01,02) distributions in the left-most four panels. Joint autocorrelations for relativistic nuclear 
collisions retain almost all correlation structure on a visualizable 2D space and provide access to 
parton fragment angular correlations with no leading-particle condition, sampling a minimum- 
bias parton distribution. Jet correlations are thus revealed with no a priori jet hypothesis, 
providing access to the low-Q 2 partons which serve as Brownian probes of the QCD medium. 

4. The p-p reference system 

The reference system for low-Q 2 partons in A-A collisions is the hard component of correlations 
in p-p collisions. The single-particle pt spectrum for p-p collisions can be decomposed into soft 
and hard components on the basis of event multiplicity dependence (Ej. Event multiplicity 



determines statistically the fraction of p-p collisions containing observable parton scattering 
(hard component). Hard components for ten multiplicity classes in the first panel of Fig. |31 
obtained by subtracting fixed soft-component spectrum model So, are plotted on transverse 
rapidity y t = ln{ (mt + Pt)/mo\. The approximately gaussian distributions on y t may be 
compared with conventional fragmentation functions plotted on logarithmic variable £ = 
\n{Ej et /pt\ \n\- Such single-particle structures motivated a study of two-particle correlations 
on (yti,yt2)- An example in Fig. |3] (second panel) reveals structures at smaller and larger y t . 




Figure 3. Single-particle hard-component distributions on transverse rapidity for ten 
multiplicity classes of 200 GeV p-p collisions; two-particle distribution on transverse rapidity and 
soft and hard components respectively of joint autocorrelations on pseudorapidity and azimuth. 

Soft and hard correlation components on yt, interpreted as longitudinal string fragments 
(smaller yt) and transverse parton fragments (larger yt), produce corresponding structures in 
joint angular autocorrelations on (?7a>0a)- in the third panel, string-fragment correlations 
for unlike-sign pairs are determined by local charge and transverse-momentum conservation 
(the sharp peak at the origin is conversion electrons). Minimum-bias parton fragments in the 
fourth panel produce classic jet correlations, with a same-side (tja < tt/2) jet cone at the 
origin and an away-side (77A > tt/2) ridge corresponding to the broad distribution of parton- 
pair centers of momentum. Similar-quality parton fragment distributions on (??Ai0a) can be 
obtained for both pts of a hadron pair down to 0.35 GeV/c (parton Q/2 ~ 1 GeV). The criteria 
for partons as Brownian probes are 1) Q 2 large enough that resulting hadron correlations are 
statistically significant and uniquely assigned to parton fragments, and 2) Q 2 small enough 
that correlations are significantly modified by local medium dynamics. In the QCD context the 
medium itself is formed from low-Q 2 partons. In the low-Q 2 regime 'partons' may not interact 
as point color charges, and complex couplings to the medium, e.g., tensor components of the 
velocity field (Hubble expansion), may be important. Non-perturbative aspects of low-Q 2 parton 
collisions should be accessible via low-p t fragment angular autocorrelations and two-particle yt 
distributions. 

5. (pt) fluctuations and prehadronic temperature/ velocity structure 

Event-wise (pt) fluctuations generally result from local event-wise changes in the shape of the 
single-particle pt spectrum, as illustrated in Fig. 0] (first panel). In each collision, a distribution 
of 'source' temperature and/or velocity on (r/, (f>) determines the local parent pt spectrum shape. 
Each hadron pt samples a spectrum shape determined by the sample location, as shown in Fig.0] 
(second panel). The local parent shape can be characterized schematically by parameter f3(rj, 4>), 
interpreted loosely as 1/T or v/c for the local pre-hadronic medium. Variation of either or both 
parameters relative to an ensemble mean results in (pt) fluctuations. 

A similar situation is encountered in studies of the cosmic microwave background (CMB) 
as shown in Fig. 0] (third panel) [T^]. The temperature distribution (3 on the unit sphere is 
represented by the microwave power density (local spectrum integral rather than mean). The 
(3(6, cj)) structure for that single event is directly observable due to large photon numbers. In 



contrast, for a single heavy ion collision as in Fig. 0] (fourth panel) the parent distribution is 
sparsely sampled by ~ 1000 final-state hadrons, and parent properties are not accessible on an 
event-wise basis. 
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Figure 4. Sources of (p t ) fluctuations, temperature/velocity distribution on binned angle 
space, WMAP microwave distribution and particle distribution from central Au-Au collision. 



Interpreting (pt) fluctuations has two aspects: 1) study equivalent two-particle number 
correlations on pt or yt, which reveal medium modification of the two-particle parton fragment 
distribution — those correlations are directly related to a distribution on (/3i,/?2) sensitive to 
in-medium parton dissipation; 2) invert the scale or bin-size dependence of (pt) fluctuations 
to obtain pt autocorrelations on {n, (ft) which reveal details of event-wise f3(n, (ft) distribution. 
We first consider properties of j3{n, (ft) as a random variable and its relation to two-particle 
correlations on pt or yt. We then employ pt autocorrelations from ijpt) fluctuations to infer 
aspects of the (3{n, (ft) distribution which depend only on separation of pairs of points on (rj, (ft). 



6. Parton Dissipation in the A-A Medium 

(pt) fluctuations can be related to a ID distribution on temperature/ velocity parameter (3 and 
corresponding two-point distribution on (^1,^2)- Each entry of those distributions corresponds to 
an event- wise pt spectrum in a single bin or pair of bins on (n,(ft). The frequency distribution on 
represents variation of the single-particle pt spectrum shape. For Gaussian-random fluctuations 
the relative variance of the j3 distribution is <t|//3q = 1/n, where n is the exponent of Levy 
distribution Aj (l+/?o(m t — mo)/n) n describing the average pt spectrum shape |19] . The shape of 
the single-particle spectrum is thus related to the event-wise temperature/velocity distribution. 
Other aspects of shape determination, such as collective radial flow, also contribute to exponent 
n. We therefore consider the two-point distribution on (/3i,/?2). 

Given the correspondence between the fluctuation distribution on (3 and the shape of the 
single-particle spectrum on p t we seek the relation between the distribution on (/3i,/32) and 
the shape of the two-particle distribution on {pa,Vt2)- The distribution on (/?i,/?2) provides 
information about the correlation structure of event-wise (3 distributions. The two-particle Levy 
distribution on (pn,Pt2), constructed as a Cartesian product of two single-particle distributions 
with Levy exponent n, represents a mixed-pair reference distribution (pairs from different but 
similar events). We can also define a two-particle object Levy distribution representing sibling 
pairs (pairs formed from single events), with exponents nj; and «a representing variances on 
sum and difference axes (/?e,/3a)- The ratio of object and reference distributions reveals a 
saddle-shaped structure whose curvatures measure temperature/velocity correlations on (n,(ft). 

Ratios of sibling to mixed pair densities for 130 GeV Au-Au collisions are shown in Fig. El 
(first two panels) plotted on variable X(pt) |20j . Those panels are dominated by a Levy saddle, 
a 2D manifestation of two-particle pt spectrum shape variation due to velocity and temperature 
fluctuations in the parent distribution. The saddle is an intermediate shape in the dissipation 
process; its curvatures reflect the correlation structure of the distribution, especially 

its covariance as discussed in |20| . The saddle curvatures on sum and difference variables, 



Figure 5. Left panels: Two-particle distributions on X(pt) (see text) for peripheral and central 
Au-Au collisions. Right panels: Corresponding number distributions on transverse rapidity. 

measured by 1/rts — 1/n and 1/n a. — 1/n, represent the variance excesses (beyond independent 
Pt sampling from a fixed parent) and covariance of temperature/velocity fluctuations for small- 
amplitude Gaussian-random fluctuations. For an equilibrated system the saddle would be 
flat (zero curvatures), and (3 fluctuations would be consistent with finite- number fluctuations: 
<t| e //3q = <?0 A /Po = The integral of correlations on {pti,Pt2), measured by the saddle- 

curvature difference l/n^ — l/n^, is equivalent to (pf) fluctuations measured in the corresponding 
detector acceptance (Hj. With increasing Au-Au centrality the curvature on the difference axis 
increases strongly, while that on the sum axis approaches zero |20| . 

More recently, we have transitioned from per-pair correlation measure f — 1 plotted on variable 
X{pt) to per-particle density ratio Ap/^yp re f plotted on transverse rapidity yt- We wish to 
follow, within a single context, the transition from parton fragment distributions in elementary 
collisions to correlations from parton dissipation in a bulk medium. Fig. (last two panels) shows 
Ap/^p re f on (yti,Ut2) for peripheral and central Au-Au collisions at 200 GeV. The logarithmic 
yt interval [1,4.5] corresponds to linear pt ~ [0.15,6] GeV/c. Peripheral collisions produce a 2D 
minimum-bias parton fragment distribution peaked at yt ~ 2.5 {pt ~ 1 GeV/c), similar to p-p 
collisions but without small-yt correlations from string fragmentation. As centrality increases the 
fragment distribution is transported to smaller yt and approaches a shape corresponding to the 
Levy saddle on X(pt) x X(pt). In this format we can study the transition with A- A centrality 
between two extreme cases: 1) in vacuo distributions of string and parton fragments and 2) 
gaussian-random variation of (3 on (77, for a nearly-equilibrated system. Parton dissipation in 
the A- A bulk medium is represented by the transition between those extremes. 

7. (p t ) fluctuations and p t autocorrelations 

The previous section describes {pt) fluctuations in terms of two-particle number densities on 
(pn,Pt2) or its logarithmic equivalent (yti,yt2), the issue being modification of the two-particle 
parton fragment distribution with changing A- A centrality. One can also express (pt) fluctuations 
in terms of two-particle pt distributions on (77, (p) which reveal different aspects of the underlying 
two-particle number distribution on vector momentum. This section describes a procedure 
to determine the correlation structure of the (3(rj,(j)) distribution as a temperature/velocity 
distribution on the prehadronic medium. 

Fluctuations in bins of a given size or scale are determined by two-particle correlations with 
characteristic lengths less than or equal to the bin scale. By measuring fluctuation magnitudes 
as a function of bin size one can recover some details of the two-particle correlation structure — 
those aspects which depend on the separation of pairs of points, not on their absolute positions. 
The relation between fluctuations and correlations is given by the integral equation 

a 2 / ^ „ A P(Pt ■ n;ke v ,le^) 

A °p t ,n\ me ni ne <t>) = 4 2^ e V e <P K mn;kl , , (1) 

k,i=i J Pref(n;ke v ,l e^) 



with kernel K mn -ki = (m — k + l/2)/m - (n — l + l/2)/n representing the 2D macrobin system, 
A<jp t . n (57],5(j)) is a variance excess and Ap(p t : n; )/^J p re f(n) is an autocorrelation density ratio. 
That equation can be inverted numerically to obtain the pt autocorrelation. 




Figure 6. Left panels: (pt) fluctuation scale dependence and corresponding p t autocorrelation 
for peripheral Au-Au collisions. Right panels: Similar distributions for mid-central collisions. 

Fig. shows fluctuation scale dependence on bin sizes (5rj,8(f)) and joint pt autocorrelations 
on difference variables (r/A, ^a) for peripheral (left panels) and mid-central (right panels) Au-Au 
collisions. Fluctuation measurements at the full STAR acceptance correspond to the single points 
at the apex of the distributions on scale. Measurements with different detectors correspond to 
different regions of those surfaces. Inversion to autocorrelations provides physical interpretation 
of fluctuation scale dependence. By inverting {pt) fluctuations, parton fragment distributions are 
visualized as temperature/velocity structures on (77, </>) complementary to number correlations on 
(yti,yt2) described in the previous section. A more comprehensive picture of parton scattering, 
dissipation and fragmentation is thereby established. 




Figure 7. pt autocorrelation from inversion for peripheral Au-Au collisions, autocorrelations 
from pair counting for minimum-bias p-p collisions and from collisions with n c h > 9. 

Pt autocorrelations can also be determined directly by pair counting. In Fig. [7| the peripheral 
Au-Au result from the previous section (first panel) is compared to the minimum-bias p-p result 
(second panel) and to p-p collisions with n c ^ > 9 (third panel). The last panel shows the 
charge-dependent (like-sign — unlike-sign pairs) pt autocorrelation for the same event class, 
reflecting charge-ordering along the jet thrust axis during parton fragmentation. This is the 
first determination of pt correlations in p-p collisions. 

8. Local velocity structure and same-side recoil 

Whether derived from pair counting or from fluctuation inversion, the resulting pt 
autocorrelations can be separated into several components. We first subtract multipoles on 
azimuth (azimuth sinusoids independent of pseudorapidity) , revealing structure associated with 
parton scattering and fragmentation. Fig. |H1 shows the resulting pt autocorrelation for 20-30% 
central Au-Au collisions at 200 GeV (first panel) and a three-component model fit to that 



distribution (second panel) including a same-side (0a < tt/2) positive peak, a same-side negative 
peak and an away-side (</>a > tt/2) positive peak. The fit is excellent, with residuals at the 
percent level. The third panel shows the result of subtracting the positive same-side model peak 
(representing parton fragments) from the data in the first panel. The shape of the negative same- 
side peak is very different from the positive peak; there is thus negligible systematic coupling in 
the fit procedure. The fourth panel shows the data distribution in the third panel plotted in a 
cylindrical format, suggesting an interpretation in terms of temperature/velocity correlations. 




Figure 8. pt autocorrelation for mid-central Au-Au collisions, model fit, data autocorrelation 
with positive same-side model peak subtracted and the same distribution in cylinder format. 

Histogram values of the pt autocorrelation effectively measure correlations (covariances) of 
blue or red shifts of local pt spectra relative to the ensemble mean spectrum at pairs of points 
separated by (t]a, (Pa)- The negative same-side peak can therefore be interpreted as a systematic 
red shift of local pt distributions adjacent to the positive fragment peak. The red shift can in 
turn be interpreted as recoil of the bulk medium in response to stopping the parton partner of the 
observed parton (positive same-side peak). This detailed picture of parton dissipation, stopping 
and fragmentation in a A-A collisions, including recoil response of the dissipative bulk medium 
suggested in the fourth panel, is accessed for the first time with joint pt autocorrelations. 

9. Reconstructing Event- wise Temperature/ velocity Structure 

We now consider the relation of pt autocorrelations to individual collision events. In Fig. El we 
repeat the WMAP CMB distribution of microwave power on the unit sphere, picturing a single 
Big Bang 'event' which has a large statistical depth and can therefore be directly observed |18j . 
Information relevant to cosmological theory is extracted as a power spectrum on polar angle 
(second panel), formally equivalent (within a Fourier transform) to an autocorrelation according 
to the Wiener-Khinchine theorem. In some studies, CMB angular autocorrelations and cross- 
correlation have been determined directly |21j . 




Figure 9. CMB power distribution, corresponding power spectrum, joint pt autocorrelation for 
mid-central Au-Au collisions and simulated temperature/velocity distribution for single collision. 

In our study of heavy ion collisions we obtain angular pt autocorrelations as in the third 
panel. Due to sparse sampling we cannot directly visualize the temperature/velocity structure 



of individual collision events as for the CMB survey. The local microwave power density of 
the CMB survey is analogous to local (pt) in a Au-Au collision. For individual collisions, and 
especially for smaller bin sizes, the event-wise mean values are not significant. However, given 
Pt autocorrelations we can simulate event-wise velocity/temperature distributions. We estimate 
the number of hard parton scatters within the STAR acceptance in a central Au-Au collision 
as 20-40, based on an analysis of p-p collisions [TS|. Combining that frequency estimate with 
shape information from the autocorrelation, and introducing some statistical variation of peak 
structure about the autocorrelation mean value, we can produce simulated events as shown in 
Fig. (fourth panel): distributions on primary angle variables (rj, whereas the autocorrelation 
is on difference variables (^A)0a)- This exercise illustrates that while Au-Au collisions are 
RHIC may be locally equilibrated prior to kinetic decoupling, they remain highly structured due 
to copious parton scattering which is not fully dissipated. Access to that structure requires pt 
rather than angular or number autocorrelations on (77, <p) to provide the full picture. 

10. Energy dependence of (p t ) fluctuations and parton scattering 

Given this close connection between parton scattering and fluctuations, the collision-energy 
dependence of (pt) fluctuations may reveal previously inaccessible parton dynamics at lower 
collision energies. In Fig. ^3 (first panel) we show the centrality dependence {y measures mean 
participant path length in nucleon diameters) of (pt) fluctuations for four RHIC energies and a 
summary (crosshatched region) of SPS fluctuation measurements at 12.6 and 17.3 GeV [22], ah 
at the full STAR acceptance (CERES measurements are extrapolated). In the second panel the 
pseudorapidity scale dependence of fluctuations at full azimuth acceptance is shown for central 
collisions at six energies. Extrapolation of CERES data in the first panel is illustrated by the 
dashed lines at the bottom of the second. Fluctuation measure Aa pt:n is related to the variance 
difference in Eq. (^Q) by A(jp t:n = 2a Pt Aa Pt:n , with ap t the single-particle variance. To good 
approximation Aa pt - n ~ <& pt , and both are per particle fluctuation measures. <£ pt was used for 
the CERES fluctuation measurements. 




Figure 10. Centrality, pseudorapidity scale and energy dependence of (pt) fluctuations for 
central collisions in the STAR acceptance and energy dependence for Hijing Monte Carlo. 

For either measure we observe a dramatic increase in {pt) fluctuations from SPS to RHIC 
energies. The centrality dependence in the first panel suggests that fluctuations for p-p and 
peripheral A-A collisions saturate near 60 GeV, whereas there is monotonic increase for the more 
central collisions. The scale dependence in the second panel illustrates how measurements with 
different detector acceptances are related. Measurements over common scale intervals should 
correspond. At RHIC energies we have demonstrated that (pt) fluctuations are dominated by 
fragments from low-Q 2 parton collisions. The energy dependence of Aa pt:n or <J> Pt is shown in 
the third panel of Fig. El plotted vs y/SNN- We observe that (pt) fluctuations vary almost 
linearly with logj-^/SAw/10.5} (solid curve in that panel), suggesting a threshold for observable 
parton scattering and fragmentation near 10 GeV. 



Fluctuation measurements based on T, pt ~ J Aa^ t . n /(n c hP^) [Hj appear to contradict 
the results described here, implying instead negligible energy dependence of (pt) fluctuations 
from SPS to RHIC. We observe that nuclear collisions at RHIC are dominated by local 
temperature/velocity structure from hard parton scattering. T, pt is a per pair measure which 
averages the local pt correlation structure dominating RHIC collisions over the entire detector 
acceptance, resulting in apparent reduction of correlations with increasing A-A centrality as 
^/^participant (p er the central limit theorem) and consequent insensitivity to contributions from 
hard scattering. We want to study separately the changes in pt production (T) and in the 
correlation structure of that produced pt (ST) prior to hadronization. H pt by construction 
estimates relative temperature fluctuations of the form 5T/T. It thus divides the structure 
problem by the production problem, greatly decreasing sensitivity to each. 

11. Summary 

We have demonstrated that low-Q 2 partons, accessed here for the first time by novel analysis 
techniques including joint autocorrelations, serve as Brownian probes of A-A collisions, being the 
softest detectable dynamical objects which experience QCD interactions as color charges. Our 
analysis of p-p correlations provides an essential reference for A-A collisions. Inversion of the 
scale dependence of (pt) fluctuations provides the first access to pt autocorrelations which reveal 
a complex parton dissipation process in A-A collisions relative to p-p collisions. We observe 
possible evidence for bulk-medium recoil in response to parton stopping. We also observe strong 
energy dependence of (pt) fluctuations, which is to be expected given the dominant role of 
scattered partons in driving those fluctuations. 
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